Steady-state entanglement of Bose-Einstein condensate and a nanomechanical 

resonator 
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We analyze the steady-state entanglement between Bose-Einstein condensate trapped inside an op- 
tical cavity with a moving end mirror (nanomechanical resonator) driven by a single mode laser. The 
quantized laser field mediates the interaction between the Bose-Einstein condensate and nanome- 
chanical resonator. In particular, we study the influence of temperature on the entanglement of the 
coupled system, and note that the steady-state entanglement is fragile with respect to temperature. 
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INTRODUCTION 



Optical nano-mcchanical systems that couple optical 
degree of freedom to the mechanical motion of a can- 
tilever have been subject of increasing investigation [1]. 
In these optomechanical systems coupling is obtained 
via radiation pressure inside a cavity [2-4], or via quan- 
tum dots [5] or ions [6]. Recently, it is made possible 
to couple mechanical resonators with the ensembles of 
atoms, where, the interaction is mediated by the field in- 
side the cavity which couples the mechanical resonators 
to the internal degrees of freedom of the atoms [7, 8], 
or to motional degrees of the freedom of the atoms [9] 
causing effects, (e.g, cooling of the mechanical resonator 
via bath of atoms [10]). In quantum meteorology, vari- 
ous targets, such as, measurement of displacement with 
larger accuracy [11] and the detection of gravitational 
waves [12] are set as milestone achievements. Engineer- 
ing entanglement in nano-mechanical systems is useful in 
quantum technologies [13]. The possibility of entangling 
the electromagnetic field with motional degree of free- 
dom of mechanical systems have been explored in var- 
ious approaches [14, 15]. In this paper, we consider a 
hybrid optomechanical system which consists of a Bose- 
Einstein condensate (BEC) trapped inside a Fabry- Perot 
cavity with a vibrating end mirror (nano-mechanical res- 
onators) driven by a single mode optical field. The intra- 
cavity field mediates the interaction of nanomechanical 
resonator with collective oscillations of the atomic den- 
sity. Hence, the motional degrees of the nano-mechanical 
resonator indirectly couples to motional degrees of the 
freedom of the Bose-Einstein condensates via optical field 
inside the cavity. Therefore, intracavity field acts as non- 
linear spring between collective atomic density and nano- 
mechanical resonator. We show that (i) the mechani- 
cal vibration of the nanomechanical resonator is entan- 
gled to the motional degree of the freedom of the BEC 
in the steady state, (ii) Furthermore, we analyze the 
steady-state entanglement as a function of temperature, 
coupling strength between BEC and field, moving mir- 
ror and field, and power of input driving laser field. In 



section 2, we model the system and explain its interac- 
tion with the environment. In section 3, we calculate the 
quantum Langevin equations for our system and solve 
the dynamics. In section 3, we quantify the steady-state 
entanglement between mechanical resonator and Bose- 
Einstein condensate. Later, in section 5, we explain the 
possibility to experimentally measure the generating en- 
tanglement. Finally, in section 6, we provide concluding 
remarks. 



2. THE MODEL 

We consider a Fabry-Perot cavity with a moving end 
mirror driven by a single-mode optical field of frequency 
up, and BEC of iV-two level atoms are trapped inside 
the Fabry-Perot cavity [16, 17]. The condensate atoms 
placed in the cavity observe one-dimensional optical lat- 
tice, formed by the oppositely propagated electromag- 
netic field inside the cavity. We consider that the atom- 
field detuning A a is very large, spontaneous emission is 
negligible, and, as a consequence we adiabatically elim- 
inate the internal excited state dynamics of the atoms. 
In addition, we also consider that the atomic densities 
are low enough that one can neglect the two-body in- 
teractions. In a weakly interacting regime, the recent 
experiment [17] suggests that only the first two symmet- 
ric momentum side modes are excited with momentum 
± 2ftk, where k is the wave number of the field. Moreover, 
we assume that the vibrating end mirror of the optical 
cavity of length L is performing harmonic oscillations 
with frequency w m along the cavity axis (x-axis). 

The Hamiltonian of the hybrid optomechanical sys- 
tem formed by the BEC, the intracavity field, and the 
vibrating end mirror of the cavity (nanomechanical res- 
onator) in the rotating frame at the laser frequency u p is 
given by [18, 19] 
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where, H a = hua.af&, Hb = fca> m b'b and H c = fiAc'c- 
i?iE(c — fiT). Here c(&) is the annihilation (creation) 
operator of the single-mode optical field, A = lo c — oj p + 
NUo/2, where U is the optical lattice barrier depth per 
photon and represents the atomic back action on the field 
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[20], and — ihF,(c — fit) shows the driving of the cav- 
ity field with amplitude E related to the laser power P 
by |E| = y/2P n/huip, where k is the decay rate of the 
photons in to the associated outgoing modes. Moreover, 
a (at) is the annihilation (creation) operator of the BEC 
mode with frequency Lo a = 4uj t , where lo t is the recoiled 
frequency. In addition, u> m and b (bt) are frequency and 
annihilation (creation) operator of the nano-mcchanical 
resonator, respectively. We assume that the frequency of 
the nanomechanical resonator is less than the free spec- 
tral range, (i.e, U) m « v/2L (v is the speed of light). 
Therefore, scattering of photons into others modes, ex- 
cept the driven mode, is neglected [19]. Here, H mc ac- 
counts for the interaction between nano-mechanical res- 
onator and intracavity field and is given by 



H„ 



-ih-^r 6*6 (b 



b t). 



(2) 



Here, the coupling strength between nano-mechanical 
resonator and light radiation pressure is defined by g mc = 
\/2(wp/L) x Q , where x D — \J h/2 m ui m is the zero-point 
motion of the mechanical mode of mass m. Moreover, 
the interaction of the BEC with intracavity mode is de- 
scribed by the Hamiltonian H ac given by 



H ac = ih — ^= 6*6 (a + a*) , 
v2 



(3) 



where, g ac = (U vN)/2 is described the strength of in- 
teraction between BEC mode and intracavity field. 

In order to describe the complete dynamics of the 
system we include the dissipation effects. In addition to 
the dynamics described by the Hamiltonian in Eq.(l), 
the system is exposed to the random noise forces due to 
quantum fluctuations of the radiation field and fluctu- 
ation of the phononic heat bath associated to the me- 
chanical resonator. We neglect the thermal effects of the 
atomic cloud and assume that the vacuum noise associ- 
ated with the cavity field is Markovian in nature with 
decay rate k and noise operator C{ n (t) of the input field 
which obeys the following correlation functions, 

(oUt)cW)) = (n c + l)5(t-t<). (4) 

Here, n c = [exp{?i ui c /Kb T} — 1] 1 is the equilibrium 
occupation number of the optical oscillator. For optical 
frequency ui c we consider Huj c /Kb T >> 1 and set n c = 0. 

The motion of the nano-mechanical resonator is 
affected due to thermal bath is Brownian and non- 
Morkovian in nature [8]. The quantum effects on me- 
chanical resonator are only observed in the limit of very 
high mechanical quality factor Q = uj m /j >> 1, the 
Brownian noise operator can be modeled as Markovian 
with the decay rate of the mechanical mode is 7. There- 
fore, the noise operator h{ n (t) can be characterized as 



where, n t h = [cxp{ftu; m / Kb T} — 1] 1 is the equilibrium 
thermal occupation number of the mechanical resonator. 



3. HEISENBERG-LANGEVIN EQUATIONS 

In order to describe the complete dynamics of the sub- 
systems involved in this problem, an adequate choice is 
to use the formalism of the quantum Langevin equations. 
Therefore, the Heisenberg-Langevin equation of motion 
for the intracavity mode, mechanical mode and bosonic 
field operator can be written as 

■ gac + 

a = — jw a a- i — = c'c, 
V2 

b = — iu m b + i c*c — 7 b + ^^27 bj n , 
v2 

_ i A + ii^(b + bt)-ii^(a + a+)-K > ) c 
v2 V2 / 

+E + V2K,c in , (6) 

where dot denotes the time derivatives and for simplicity 
we omit the hat symbol from the operators. These are the 
nonlinear quantum Langevin equations and dynamics is 
complicated. In the following we linearized the operators 
around the steady state values, a = (a) ss + 9a, b = (b) ss + 
9b, c = (c) ss + 9c. Here, we assume that the fluctuation 
operators 9a, 9b and 9c have zero mean. The steady 
state value of the intracavity mode is (c) = E/(k + zA), 
where the total effective detuning is 
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\ / SS 
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(7) 



For the sake of simplicity we assume that the field 
is real positive and this can be achieved by adjust- 
ing the phase of the laser field. Similarly, the steady 
state value of the BEC and mechanical resonator 
modes are (a) ss = [-g ac /\/2w a ] (c+c)^ and (b) sa = 
[i g mc /V2(7 + i w m )] (etc ) sg respectively. 

We linearize the Langevin equations of motion given 
in Eq.(6), and assume that pump field is intense and keep 
terms only up to first order in the fluctuation operators. 
We rewrite each Hcisenberg operator in Eq.(6) as a sum 
of steady state value and fluctuation operator with zero 
mean value. Therefore, the linear set of equations are, 

9a = — iw a 9a — *— r- (9c + 9ct) , 

9b = — (7 + iuj nl )db + *— ^— (9c + 9ct) - 

96 = -( K + i A)9c + z^p (9b + 9bt) 
.G a 



27 bm 



"*~<p ip & + ^ a ^) + v2Kc i: 



(8) 



(5) 



The linearized quantum Langevin equations show that 
the fluctuations of mechanical resonator and BEC are 



3 



now coupled to the cavity field quadrature fluctuations 
by the effective couplings G mc = \/2gmc (c) ss and G ac = 

\/2gac ( C ) SS ' wn i cn can be made very large by increasing 
the amplitude (c) of the intracavity field. Linearized 
quantum Langevin equations (8) and their correspond- 
ing Hermitian conjugate form a system of six first-order 
coupled operator equations, for which the Ruth-Hurwitz 
criteria implies that the system will be stable only if the 
following stability condition is satisfied 



(A 



K 2 ) W m W a 



AG 



2 

mc 



Gar 



> 0. 



(9) 



In the following we consider only red-detuning regime 
(A > 0) and from now on we assume that the 
above stability condition is satisfied. The quantum 
Langevin equations (8) are linear in creation and annihi- 
lation operators. We transform to the quadratures, i.e, 
9q m = (db + 9bt)/V2, dp m = (db - dtf)/iV2, 8q a = 
(9a + <9at)/V2, <9 Pa = (da, - de!)/iy/2, dq c = (dc + 
dJ)/V2, dp c = (dc - dc^)/iV2. The system of lin- 
earized equations of motion can be written in com- 
pact matrix form as R(f) = MR(i) + F(t), where, 
R = (dq m , dp m , dq a , dp a , dq c , 9p c )* is the vector of 
the quadrature fluctuations, and in superscript t de- 
scribes transpose of matrix. Furthermore, F is the vec- 
tor corresponding to noises, whereas M is drift ma- 
trix. Since the quantum noises are white in nature 
and the dynamics is linearized, hence the state of the 
system will be a zero mean Gaussian state, and there- 
fore completely determined by the covariance matrix 
Vy = ((Ri Rj + Rj R) - 2R Rj). In order to find the 
steady state covariance matrix V, we solve the linearized 
quantum Langevin equations as did in [15]. In steady 
state the covariance matrix fulfills the Lyapunove equa- 
tion 



MV + VM' 



(10) 



Eq.(10) is the linear matrix equation and can straight 
forwardly be solved. However, the general exact expres- 
sion is too cumbersome and not reported here. One can 
extract all the information about the steady state of the 
system from the correlation matrix. 



4. STEADY STATE ENTANGLEMENT 

We compute the entanglement between mechanical 
mode and the BEC mode in steady-state by tracing out 
the cavity mode. The steady state entanglement is deter- 
mined by computing the logarithmic negativity En from 
the corresponding covariance matrix V. We consider an 
example where the length of the cavity L = 1 mm and 
laser input power P = 50 mW [15]. The effective mass 
and resonance frequency of the mechanical resonator are 
m = 4ng and w m = 1 MHz. The decay rate of the me- 
chanical resonator is 7 = 2ir x 100 Hz and cavity finesse 
F = 1 x 10 4 . For these values of the parameters, the cou- 
pling g mc between mechanical resonator is in the order of 





gmc< H '> 



(b) 



FIG. 1: (Color online) (a) Plot of logarithmic negativity 
En as a function of coupling of the intracavity mode with 
nano-mechanical resonator g mc and BEC mode g ac for A = 
2ir x 2 MHz. The optical cavity length L = 1 mm, driven by 
a laser with wavelength A = 1000 nm, power P = 50 mW and 
u> m ~ u) & . The mechanical mirror has a frequency uj m /2n = 
1 MHz, and damping rate 7 = 2n x 100 Hz, its temperature is 
T = 10 fjK and cavity finesse is F = 1 x 10 4 . (b) En versus 
g mc and g a c for A = 2tt x 3 MHz. 



kHz. Moreover, the interaction of the optical field with 
Bose-Einstein condensate is kept small, so that Bogoli- 
ubove mode expansion becomes possible. 

Measurement of the entanglement between the me- 
chanical resonator and Bose-Einstein condensate, re- 
quires as to compute En, which is obtained by tracing 
out the cavity mode, i. e, removing the rows and columns 
of V which correspond to the cavity mode. The reduce 
state is still Gaussian and fully characterized by 4 x 4 ma- 
trix V. In order to measure the entanglement between 
mechanical mirror and intracavity filed, we consider the 
Logarithmic negativity En- In the case of continuous 
variable (CV), En can be defined as [21] 



En = max [0, — In 2 , 



(11) 



where, v_ = 2- 1 / 2 {£(V) - [£(V') 2 - 4dct V'] 1 / 2 } 1 / 2 , 
with XX v ') = dctx + dctY - 2dctZ, is the smallest 
symplectic eigenvalue. However, the second eigenvalue 
u+ = 2- 1 / 2 {£(V) + E(V') 2 - 4det V'] 1 / 2 } 1 / 2 » 1/2 
at any value of the parameters. Therefore, it has no effect 
on the non-separability of the state [22]. Moreover, the 
correlation matrix V' in 2 x 2 block form can be written 



as 



V'= 



X z 
7} Y 



It is clear from Eq.(ff) that the En is the decreasing 
function of , f_ and it quantifies how much two Gaussian 
states are entangled. The Gaussian state gets entangled 
only if V- < 1/2, and it is Simon's necessary and suffi- 
cient entanglement non-positive partial transpose crite- 
rion of the Gaussian states [22], and this condition can 
also be written as 4det V < £(V) - 1/4. 

In Fig. 1 we show that the entanglement between 
nano-mechanical resonator and Bose-Einstein conden- 
sate quantify by Logarithmic negativity En as a func- 
tion of the coupling of the mechanical resonator and 
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Bosc-Einstcin condensate with field in the cavity. In 
Fig. 1(a), we have En — 0.15 for g mc = g ac ~ 300 Hz 
and A = 2tt x 2 MHz. In Fig. 1(b), we find that the 
E N ~ 0.2 for g mc = g ac ~ 300 Hz and A = 2tt x 3 MHz. 




50 T,nK) 



FIG. 2: (Color online) (a) Logarithmic negativity En as a 
function of power of the input laser field and normalized de- 
tuning A/u m . (b) En as a function of coupling strength be- 
tween mechanical mirror and intracavity field, i.e, g mc and 
normalized detuning A/ui m for a constant value of the cou- 
pling strength between BEC and intracavity field, i.e, g ac = 
100 Hz. (c) En as a function of A/cj m for m = 4ng(black 
solid), m = 5ng(red dashed) and g ac = 100 Hz. (d) En as 
a function of temperature T for A ~ 4cj m , m = 4ng (black 
solid), m = 5 ng (red dashed) and g ac = 100 Hz. All the others 
parameters are the same as in Fig. 1. 

In Fig. 2, we also study the robustness of the steady 
state entanglement of the nano-mechanical resonator to 
the Bosc-Einstcin condensate with respect to several pa- 
rameters in system. Fig. 2(a) shows the dependence of 
En on the input laser power P and normalized detun- 
ing A/oi m for g ac = 100 Hz. Fig. 2(b) shows the depen- 
dence of En on optomechanical coupling rate g mc and 
normalized detuning A/uj m for g ac = 100 Hz. Moreover, 
in Fig. 2(c) we plot the En as function of normalized de- 
tuning A/w m for two different values of the mass m of 
the mechanical resonator, i.e, black solid and red dashed 
lines for m = 4 ng and m = 5 ng respectively. One can 
easily observe that En exists only within a finite interval 
of values of the A around A ~ 4.5 w m . Now, we inves- 
tigate the effects of temperature on the entanglement. 
We consider En of the evolved system at a fixed value of 
detuning. At high temperature, thermal fluctuation al- 
ways suppress the entanglement of the coupling systems. 
Therefore, the logarithmic negativity is the decreasing 
function of temperature. In Fig. 2(d), we plot the loga- 
rithmic negativity as a function of temperature at a fixed 



value of the detuning A = 2tt x 4 MHz. The black solid 
line refers to mass m = 4ng and red dashed line refers 
to mass m = 5ng. Fig. 2(d) shows that En monotoni- 
cally decreases with temperature. Moreover, it is noted 
that En is also very sensitive to the mass of mechanical 
resonator and quickly decay for large value of the mass. 

Latter, we find the effective coupling between me- 
chanical resonator and Bosc-Einstcin condensate. We 
consider the regime in which the dynamics of the cav- 
ity mode remains unperturbed due to the motion of the 
mechanical mode and BEC mode, and slowly mediates 
the interaction between the two. Therefore, we chose the 
regime where the cavity mode can be eliminated adiabat- 
ically. The Hamiltonian corresponding to the linearized 
quantum Langevin equations for the fluctuations opera- 
tor <9q m , <9p m , <9q a , <9p a ,<9q c and <9p c is given by 



2 

-G„ 



I = ^ (dpi + 

■ dq c <9q m + 



G ac <9q c <9q a 



(12) 



The corresponding parameteric regime for the fast cavity 
mode dynamics is, A >> G mc ,G a c or k >> G mc ,G a c- 
Due to the second condition, cavity-mediated coherent 
dynamics is destroyed. Therefore, we only consider the 
regime where A takes large values, and the fluctuation 
quadratures of the cavity mode adiabatically follow the 
dynamics of positions fluctuations of the mechanical res- 
onator and BEC mode. We further assume that both 
the mechanical resonator and BEC modes to be on res- 
onance, i.e, u; a = u> m = bj. On eliminating the photon 
degree of freedom adiabatically we get the following ef- 
fective Hamiltonian 



Heft 



tin 

2 



dpi 



<9q m <9q a , 



(13) 



where, wi = 4 G^ c A/(k 2 + A 2 ), lo 2 = 4 G 2 C A/(k 2 + A 2 ) 
and G ma = -8 G ac G mc A/(« 2 + A 2 ). It is noted that the 
effective interaction G ma between mechanical resonator 
and Bose-Einstein condensate is increased as the coupling 
of mechanical oscillator and BEC is increased with the 
intracavity field. The effective interaction between a me- 
chanical resonator and Bosc-Einstein condensate can be 
described via the Hamiltonian "H ma , from Eq.(13) 

n ma = G ma (9a9b t +9a t 0b)/2 + G ma (0a0b + aa t 0b t )/2, 

(14) 

where G ma is the effective coupling strength between me- 
chanical mode and BEC. This Hamiltonian is analogous 
to the interaction of two optical fields in a nonlinear 
medium generating the 2-modes squeezed state [24] . The 
first term in Eq.(14) corresponds to an energy exchange 
take place between mechanical and atomic modes. The 
second term in Eq.(14) accounts for the down-conversion 
interaction, describes the creation and annihilation of the 
atomic and mechanical modes phonons in pairs, which 
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corresponds to entangling the atomic and mechanical 
modes. For G ma ~ cj m or greater, the mechanical res- 
onator and Bosc-Einstein condensate are entangled, and 
the parameters are chosen such that this condition is ful- 
filled. 



5. EXPERIMENTAL DETECTION 

For experimental realization of the generated entangle- 
ment, one has to measure several quadrature correlations 
[23] as have been experimentally measured for the entan- 
glement of the two optical modes [25]. However, in our 
case we consider another Fabry Perot cavity adjacent to 
first one and is driven by a weak laser field, a scheme of 
this kind has been discussed in [15]. In addition, we as- 
sume that the movable mirror is perfect reflector at both 
sides so there is entanglement between the optical modes 
of the two cavities. The equation of motion of the an- 
nihilation operator Ci of the optical mode of the second 
cavity similar to the linearized version of Eq.(8) is 



dc i 



G 



(ki + iA 1 )da + (db + 0bt) 



-i acl (<9a + dsJ) 



^inl 



(15) 



where Kx, Ai and C; ni are the cavity decay rate, effective 
detuning and the input noise of the second cavity mode 
respectively. In addition, G mc i and G ac i arc the effective 
coupling rates of the second cavity mode to the mechani- 
cal resonator and Bose-Einstein condensate respectively. 
Moreover, we assume (c) ss >> (c) ssl and Ai = u) m » k, 
Gmci, Gael- Therefore, in the rotating frame at Ai = uj m , 
Eq.(15) for slow variables do(t) = do(t)exp {iuj m t) can be 
rewritten as 




Gmci 9b - G ac i<9a 



(16) 



Here, we assume w m = LJ a = ui and neglect the fast oscil- 
lating terms at frequency 2w. On eliminating the cavity 
mode adiabatically, Eq.(16) can be written as 



dc i 



2«i 



G mc i<% - G ac i<9a 



Cinl- 



(17) 



According to the stranded input-output relation [26], 
Couti = \/2 ki dci — Cj„i, the output field is given by 



2kx 



G mc i db - G a ci 95, 



+ Cinl. 



(18) 



From Eq.(18) we can measure the G mc i<96— G ac i<9a by ho- 
modyning the output of the second cavity mode. There- 
fore, all the entries of the correlation matrix can be deter- 
mined by measuring the correlation between the output 
of the two cavities, hence from the correlation matrix 
logarithmic negativity can be calculated numerically. 
6. CONCLUSION 

In conclusion, we discuss quantum correlation in a sys- 
tem which consists of Bose-Einstein condensates trapped 
inside a Fabry-Perot cavity with a moving end mir- 
ror driven by single mode optical field. We describe a 
scheme to generate the steady state entanglement of the 
motional degree of freedom of the nano-mechanical res- 
onator (end mirror of the cavity) and collective oscillation 
of the atomic density of the Bose-Einstein condensate. 
Moreover, it is observed that the entanglement generated 
between Bosc-Einstcin condensate and nano-mcchanical 
resonator is very sensitive with respect to temperature 
and persist only upto 20 yuK and 100 ^iK for mass 4ng 
and 5 ng, respectively. The generated entanglement is 
measured by considering the second optical cavity driven 
by a weak laser field. 
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